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Ultra  Sensitive  Mass  Sensors  based  on  Suspended  Carbon  Nanotubes 


Final  research  report  of  Yaish  group 


General- 

This  project  aims  to  utilize  suspended  Carbon  Nanotube  (NTs)  as  ultra  sensitive 
mass  sensors.  For  that  purpose  few  steps  had  to  be  addressed.  The  most  important  and 
difficult  mission  was  to  fabricate  suspended  resonators  made  out  of  NTs  which 
behave  as  field  effect  transistors  (FETs)  and  posses  noiseless  and  reproducible 
electrical  characteristic.  Without  these  requirements  this  project  can't  proceed. 
Therefore,  first  we  developed  and  construct  our  own  chemical  vapor  deposition 
(CVD)  system.  We  design  and  build  the  plumbing  for  the  machine  and  install  it  in  our 
facilities.  We  succeed  in  operating  it  and  receive  high  yield  of  single  wall  and 
multiwall  NTs.  We  fabricate  devices  using  these  tubes  and  received  FET  behavior  as 
expected.  A  huge  effort  was  addressed  in  order  to  obtain  straight  and  noiseless  device. 
For  that  purpose  we  develop  novel  growth  process,  and  study  the  various  parameters 
that  affect  the  growth  outcome.  Two  kinds  of  lithography  were  implemented  and 
different  metal  and  annealing  conditions  were  tested.  At  the  end  of  this  tedious 
process  we  succeed  to  obtain  long,  straight,  and  noiseless  FET  devices  based  on  NTs. 
The  next  step  was  to  design  and  build  the  experimental  setup  for  the  vibrational 
measurements.  We  designed  and  built  three  independent  setups  for  these 
measurements  and  validate  their  operational  capabilities.  The  first  setup  include  high 
frequency  room  temperature  vacuum  chamber.  The  second  consists  of  variable 
temperature  vacuum  insert,  and  the  last  is  equipped  with  strong  magnetic  field.  All 
three  setups  have  been  tested  for  operational  at  high  frequencies  up  to  -lGHertz. 
Equipped  with  experimental  setup,  the  suspended  resonators  were  fabricated.  We 
designed  and  implemented  three  different  paths  for  fabrication  of  suspended  NTs,  and 
all  three  of  them  exhibit  promising  results.  The  first  method  was  based  on  suspending 
the  NT  up  to  its  electrical  contacts.  In  the  second  a  trench  was  made  on  top  of  the  NT 
and  only  portion  of  it  was  suspended.  The  last  method  was  based  on  growing  the  NT 
above  predefined  trench.  The  best  results  that  we  obtain  were  based  on  the  second 
method.  However,  higher  quality  factors  are  expected  to  be  using  the  third  method. 
Currently  we  pursue  our  research  in  this  direction  as  well.  For  the  second  method  we 
succeed  to  obtain  the  first  few  vibrational  modes  and  receive  wide  range  of  mode 
tunability  using  external  electric  fields.  We  developed  theoretical  model  for  the 
experimental  data  with  an  excellent  agreement.  In  addition,  theoretical  simulation  of 
NT  resonators  was  developed  and  several  interesting  prediction  were  found. 

In  addition  to  NTs,  we  developed  novel  technique  to  fabricate  suspended  Silicon 
Nanowires  (SiNWs)  that  operate  as  FETs  as  well.  This  new  approach  is  very 
important  since  it  enables  us  to  implement  the  same  detection  scheme  for  NT 
resonators  to  SiNWs  as  well.  The  chemistry  of  silicon  surface  is  very  reach  and  well 
established  thus  SiNW  resonators  offer  us  variety  of  optional  antibody-antigen  pairs 
for  sensitive  mass  detection.  Additionally,  we  studied  how  to  attach  various  molecules 
onto  the  NT’s  surface.  This  capability  is  an  important  ingredient  of  this  research 
proposal  since  enables  us  to  attach  specific  molecules  and  thus  obtain  selective  and 
sensitive  mass  sensor.  There  could  be  various  means  of  attachment.  The  attachment 
can  be  very  rigid  and  strong,  as  covalent  bonding,  but  can  be  also  looser,  like  Van  der 
Waals  interaction.  We  studied  both  regimes,  since  in  advance  it  is  not  clear  which 
type  of  connection  is  preferable  for  our  application.  Van  der  Waals  interaction  was 


studied  with  deposition  of  the  molecules  over  the  final  device.  However,  covalent 
bonding  requires  intensive  research  even  before  the  final  device  is  available. 
Therefore  we  decided  to  use  amine  chemistry  in  order  to  attach  nanocrystal  or  small 
molecules  onto  the  1-dimesional  wires.  In  practice  we  utilized  this  chemical  reaction 
to  attach  gold  colloid  to  the  NTs. 

This  final  report  is  organized  as  follows: 

•  Carbon  nanotube  growth 

•  FET  fabrication 

•  Experimental  setup 

•  Vibrational  modes 

•  Theoretical  simulation 

•  Silicon  nanowires  resonators 

•  Chemical  modifications 

Carbon  nanotube  growth: 

In  order  to  measure  the  vibration  modes  of  suspended  CNT  we  first  had  to  fine- 
tune  the  growth  process  to  achieve  reasonable  yield  of  SWNT  longer  than  20um  and 
fairly  straight,  since  curves  usually  mean  defects  in  the  CNT  structure  which  tamper 
conductivity.  Different  recipes,  catalysts,  and  substrate  treatments  were  examined. 

A  summery  of  major  modification  and  their  results  is  listed  below: 


(Flying  Catalyst) 

Recipe 

1 .  Heat  the  oven  in  backward  position  to  900°C. 

2.  Growth  at  900°C  with  0.52  SLM  CH4,  0.7  SLM.  Measure  10  minutes  when  the  temperature 
reaches  820°C. 

Remarks:  No  gas  guide _ 


Recipe 

1 .  Heat  the  oven  in  backward  position  to  900°C. 

2.  Flow  0.52  SLM  CH4,  0.7  SLM  H2.  Bring  the  furnace  to  forward  position.  Measure  10  minutes 
when  the  temperature  reaches  820°C. 

3.  Pump  the  chamber  and  cool  down  with  0.1  SLM  Ar. 

Remarks:  No  gas  guide _ 


A4  (~2nm  Fe  +  Catalyst) 


H 


A2  (~2nm  Fe  +  Fe  NP  low 
cone.) 


Recipe 

1 .  Heat  the  oven  in  backward  position  to  900°C. 

2.  Flow  2  SLM  CH4,  2  SLM  H2.  Bring  the  furnace  to  forward  position.  Measure  10  minutes  when 
the  temperature  reaches  820°C. 

3.  Pump  the  chamber  and  cool  down  with  0.1  SLM  Ar. 

Remarks:  No  gas  guide 


Heat  to  750°C  in  vacuum. 

Reduction  at  to  750°C  with  0.5  SLM  H2  and  0.5  SLM  Ar. 

Bring  furnace  backward  and  cool  down  the  chamber  to  350°C. 

Heat  the  furnace  to  900°C. 

Flow  2  SLM  CH4,  2  SLM  H2.  Bring  the  furnace  to  forward  position.  Measure  10  minutes  when 
the  temperature  reaches  820°C. 

Remarks:  No  gas  guide 


A4  (Fe  evaporation  +  Catalyst) 


A5  (Catalyst) 


Recipe 

1 .  Heat  the  oven  in  backward  position  to  900°C. 

2.  Flow  2  SLM  CH4,  2  SLM  H2.  Bring  the  furnace  to  forward  position.  Measure  10  minutes  when 
the  temperature  reaches  820°C. 

3.  Pump  the  chamber  and  cool  down  with  0.1  SLM  Ar. 

Remarks:  No  gas  guide _ _ _ 

A1  (~2nm  Fe)  I  A2  (~2nm  Fe  +  catalyst)  I  B1  (Ni  evaporation) 


n  xHj.Gfe  SE(U|  5/158  ViH 


Recipe 

1 .  Heat  the  oven  in  backward  position  to  900°C. 

2.  Flow  1.2  SLM  H2 

3.  Flow  0.52  SLM  CH4,  0.7  SLM  H2.  Bring  the  furnace  to  forward  position.  Measure  10  minutes 
when  the  temperature  reaches  820°C. 

4.  Pump  the  chamber  and  cool  down  with  0.1  SLM  Ar. 

Remarks:  No  gas  guide _ _ _ 

A3  (  ~2nm  Fe  +  Catalyst)  I  A4  (  Alumina  +  Fe  NP)  I  A2 


Recipe 

1 .  Heat  the  oven  in  backward  position  to  900°C. 

2.  Flow  0.52  SLM  CH4,  0.7  SLM  H2.  Bring  the  furnace  to  forward  position.  Measure  10  minutes 
when  the  temperature  reaches  820°C. 

3.  _ Pump  the  chamber  and  cool  down  with  0.1  SLM  Ar. _ 


Recipe 

1. 

2. 

3. 


Heat  the  oven  in  backward  position  to  900°C. 

Flow  0.52  SLM  CH4,  0.7  SLM  H2. 

Bring  the  furnace  to  forward  position.  Retract  the  furnace  at  500°C. 


After  huge  amount  of  effort  we  have  managed  to  repeatedly  grow  straight  and  long 
individual  SWNT.  Devices  fabricated  from  these  NT  exhibited  very  low  contact 
resistance  (around  50  kOhm)  and  much  better  signal  to  noise  ratio  than  devices 
fabricated  from  previously  grown  NTs. 

Field  Effect  Transistor  Fabrication: 

After  growth  optimization  electrical  contact  were  fabricated.  Here,  two  methods 
have  been  used.  The  first  was  based  on  e-beam  lithography,  where  each  NT  was 
addressed  separately.  The  main  advantage  of  this  method  is  that  each  device  is  built 
from  a  single  tube,  and  its  length  and  diameter  are  known  before  completing  the 
device.  However,  this  method  is  very  tedious,  and  difficult  to  utilize  in  order  to  study 
which  for  metal  contact  is  better  and  what  annealing  conditions  one  should  apply.  The 
second  method  is  based  on  optical  lithography.  We  designed  a  clever  mask  in  which 
after  catalyst  size  and  recipe  optimization,  statistically  single  NT  is  bridging  between 
two  nearby  electrodes  to  form  electrical  contact.  This  procedure  lack  the  knowledge 
how  many  NTs  are  connected  between  two  electrodes,  as  well  as  their  lengths  and 
diameters.  However,  this  information  can  be  achieved  after  device  fabrication,  and  the 
parallel  process  enable  us  to  optimize  growth  recipes,  metallic  electrodes,  and 
annealing  conditions. 

Below  we  present  few  electrical  characteristics  for  NT  FETs  that  were  fabricated 
in  our  lab.  In  the  first  group,  devices  were  made  by  e-beam  lithography  and  in  the  last 
by  photolithography. 

SWNT  FET-transistors  made  of  semi-conducting  cut's  (left)  and  small  band  gap  cut's  (right)  grown 
before  optimal  recipe  with  Cr/Au  contacts. 


SWNT  FET-transistors  made  of  semi-conducting  cut's  (left)  and  small  band  gap  cut's  (right)  grown 
with  optimal  recipe  with  Cr/Au  contacts. 


For  the  contact  electrodes  we  have  tries  both  Cr/Au  (2nm/l  lOnm)  and  Pd  (30nm). 
Pd  exhibited  lower  contact  resistivity  but  devices  were  noisier  (see  below).  An 
annealing  step  may  solve  this  problem. 

SWNT  FET-transistors  made  of  semi-conducting  cut's  (right)  and  small  band  gap  cut's  (left)  grown 
with  Pd  contacts 


These  devices  were  made  by  e-beam  lithography.  The  one  below  utilized 
photolithography  procedure.  In  all  of  them  Pd  contact  were  deposited. 


Since  the  yield  per  process  of  the  photolithography  path  is  much  higher,  it  is  easier 
to  find  long  and  noiseless  NTs  with  compare  to  the  e-beam  method. 

Experimental  setup: 

Two  different  measurement  chambers  were  built  for  this  measurement.  The  first  is 
a  Vacuum  chamber  with  an  option  to  introduce  different  gases.  In  the  chamber  sits  a 
PCB  which  is  connected  to  external  co-ax  cables  via  28  SMB  connectors.  2  different 
PCB  were  fabricated: 

The  first  includes  a  28  pin  socket  and  transmission  line  designed  for  high 
frequency. 

Half  of  the  lines  contain  a  50ohm  resistor  in  line  with  a  InF  capacitor  which  runs 
to  ground  just  before  the  connection  to  the  socket  to  avoid  reflection  of  incoming 
wave  due  to  the  high  impedance  of  the  device.  The  PCB  was  tested  to  give  reasonable 
signal  up  to  frequencies  of  500  MHz.  A  second  PCB  is  identical  to  first  one  apart 
from  the  fact  that  it  contains  no  socket.  The  sample  is  attached  to  the  center  of  the 
PCB  using  a  silver  paint  and  the  pads  are  directly  bonded  to  the  near-by  PCB  strip 
lines.  This  enables  higher  frequency  measurements  up  to  800-900  MHz. 

Picture  of  the  Vacuum  Chamber,  Inner  low  noise  amplifier  and  low  pass  filter 


A  second  apparatus  was  built  for  cryogenic  temperatures  down  to  4K.  It  is  a  dip¬ 
stick  that  can  be  cooled  either  by  exchange  gas  or  in  high- vacuum  using  a  cold- finger. 
In  the  bottom  of  dip-stick  there  is  a  socket  identical  to  the  one  situated  on  the  PCB  of 
the  vacuum  chamber  so  samples  can  be  interchanged.  The  socket  is  connected  to  10 
stainless-steel  flexible  coaxial  cables,  and  4  high  frequency  semi-rigid  stainless  steel 
cables.  The  flexible  cables  were  tested  to  transmit  signals  up  to  60  MHz  and  the  semi¬ 
rigid  up  to  500MHz.  The  dip-stick  temperature  is  controlled  using  an  internal  heater 
and  a  calibrated  silicon  diode  for  exact  temperature  measurements.  The  sample  can  be 
heated  from  4K  to  340K  while  in  Dewar  as  was  verified  by  a  variator  bonded  to 
socket. 


The  measurement  apparatus  is  as  follows: 


NT  device 


An  RF  signal  with  frequency  oj()  if  fed  into  the  source  contact.  A  second  RF  signal 
with  a  frequency  a>0  -  A  <x>  is  fed  into  the  gate  contact.  The  mixing  term  as  explained 


below  is  proportional  to  the  derivative  of  the  conductance  (G)  with  respect  to  the  gate 
voltage.  The  derivative  of  the  conductance  (G)  with  respect  to  the  gate  voltage  can  be 
compared  with  the  numerical  derivative  of  current  versus  gate  voltage  DC 
measurement.  The  gate  capacitance  (Cg)  is  a  function  of  the  NT's  height  above 
equilibrium  (  dz  ). 


G  =  G(q) 
q  =  Cs{z)Vg 

ISD  =  VSDG  =  (VSDDC  +  VSD) 
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In  order  to  test  the  system  we  first  made  measurements  with  devices  which  are  not 
suspended,  thus  the  first  term  in  the  brackets  of  the  last  equation  disappears. 
Comparing  the  numerical  derivative  and  the  mixed  signal  revealed  that  the  signal  is 
too  noisy  to  measure  the  small  signal  from  the  mechanical  vibrations. 


Measured  Mixing  signal  taken  at  HF=10Mhz,  IF=lkHz,  Vgate(rms)=l  3mV,  Vsd(rms)=10mv. 


To  solve  the  issue  a  low  noise  amplifier  for  the  range  of  the  intermediate  frequency 
(300Hz-8MHz)  was  designed  and  built. 


0P1.0P2  AD797 


But  the  high  frequency  signal  entering  the  amplifier  saturated  it  and  masked  the 
signal.  A  lOOdB  Decade  LPF  just  before  the  amplifier  solved  the  issue. 
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The  low-noise  amplifier  scheme  (left)  and  its  frequency  response  (right). 


As  can  be  seen  the  signal  quality  improved  significantly. 


Measured  Mixing  signal  taken  at  HF=10Mhz,  IF=lkHz,  Vgate(rms)=8mV,  Vsd(rms)=6mv. 

Suspended  Device: 

Once  good  devices  were  measured  a  second  e-beam  lithography  stage  defined 
trenches  between  the  NT's  electrodes.  A  wet  etch  and  critical  point  drying  process 
was  utilized  to  remove  300nm  and  leave  doubly  clamped  SWNT  suspended  above  a 
trench.  Alternatively  pre-growth  trench  was  defined  and  NT  grows  above  the  trench. 
Below  are  images  of  two  suspended  devices  with  these  two  different  procedures. 
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Vibrational  modes: 

All  the  steps  that  we  took  in  this  research  boil  to  the  success  of  measuring  the 
vibrational  modes  of  suspended  NTs.  We  succeed  to  measure  few  resonance 
frequencies  of  2  micron  length  and  1.4nm  diameter  suspended  NT.  The  electrical 
characteristics  of  this  device  were  very  quite,  and  posses  small  band  gap 
semiconductor  NT  behavior. 
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V[V] 

When  introduce  into  the  vacuum  chamber  and  perform  our  mixing  experiment  the 
following  graph  was  obtained.  The  x-axis  represents  the  dc  gate  voltage  that  may  tune 
the  tube  and  hence  modify  its  resonance  frequencies,  and  the  y-axis  represents  the 
resonance  frequencies.  The  different  colors  present  the  tube  conductance. 
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Vg[V] 


One  can  notify  the  dim  contour  that  have  parabolic  shape  and  extend  from  low  to 
high  frequency  as  the  gate  voltage  increases.  A  cut  along  the  x-axis  is  plotted  below. 
The  arrows  point  on  the  resonance  frequencies  of  this  tube. 


In  order  to  enhance 
voltage  is  plotted  below. 


the  visibility  of  these  lines,  derivative  with  respect  to  gate 
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Vg[V] 


Now,  one  can  see  very  easily  the  dispersion  curves  of  the  vibrational  frequencies. 
In  order  to  develop  theoretical  model  for  these  resonances  one  should  address  very 
carefully  all  the  factors  that  affect  these  modes.  More  specifically,  the  current  along 
the  NT  is 

I  =  G  Vsd 
r)C 

61  =  —6q  ■  6Vad(t) 
dq 

Sq  =  CgSVg(t)  +  VgSCg 

6Vsd(t)  =  6Vsd  •  cos (ujit  +  tpi) 

6Vg(t )  =  6Vg  ■  cos(u>2t  +  <fi2 ) 

F)C 

tC.it)  =  -~-6z(t) 

All  together  we  get 
1  d  C 

iI  =  C"dV'  {C’6V^  +  V9SC9(t))  ■  SV,d(t) 

i  fin 

=  CrgWg'  (CWW  +  V9~dfSZ{t))  '  Wsd{t) 

The  current  expression  contains  5z(t).  In  order  to  know  z(t)  we  need  to  calculate 
the  force.  After,  we  can  solve  the  differential  equation  for  forced  and  damped 
harmonic  oscillators.  The  force  that  acts  on  the  NT  is: 

1  riC 

F=2  -Qf(V9DC  +  SV9(t)~5Vad(t)Y 

Since  the  DC  term  is  much  bigger  than  the  AC  term  we  can  approximate  it  to 

l  / )( ‘ 

F = m((v°DC)2 + 2v*DC{sv^  -  sv°dm 

Now  we  can  solve  the  harmonic  equation: 


dC 

m  •  z  +  T  •  z  +  k  ■  z  =  —  -  VgDC  •  (<5V^  •  cos(lj2 t  +  <£>2)  —  •  cos(u?i t  +  <^1)) 

For  conventional  geometry  where  the  NT  is  suspended  above  planar  gate  the 
capacitance  is  given  by 

_  27T£q 

9  ln(2 z/r 

dCy  1 

dz  9  z\n(2z/r) 

When  solving  the  harmonic  equation  one  finds  8z(t)  =  A  •  cosfwt  +  qp)  and  the 
current  is 

dG  1 

SI  -  -=T-(5Vgcos(uj2t+(p2)-VqDC  -  Acos(cot+y?))-8Vsdcos(uJi^+y>i) 

oVg  zln(zz/rj 

In  the  lockin  this  term  is  multiply  by  cos((g)1  -ro2)t+cp3)  and  we  perform  an 
integral  over  2tt/ (co  1  -  co2).  Lets  find  the  solution  for  5z(t).  We  ignore  for  the 
moment  that  the  force  term  has  two  frequencies  and  assume  we  need  to  find  solution 
for  the  following  equation 

mz  +  Tz  +  kz  —  fa  ■  cos(ut  +  p) 

The  solution  is 


z{t)  = 


fo 


v/m2(cUo-w2)2  +  (iy)2 


,  m(o?n  —  a/2)  ,  . 

{  — .  cos  (uft  T  <p)  T- 

^  /  0/  o  o\o  /T~I  \o  v  r  ' 


Tlo 


s/rrfifuj q  -  tv2)2  +  (Tlj)2 

_ fo _ 

Vm2(a;2-u;2)2  +  (ra;)2 


—  u>2)2  +  (Tw)2 

cos(o>t  +  p  —  p(u)) 


sin(u;f  +  9)} 

(2) 


where  we  use  the  relation  cos(a  -  b)  =  cos(a)  cos(b)  +  sin(a)  sin(b) 
and 


cos(<^(o;))  = 

sin(ip(u))  = 
Also  we  have 
tan(</?(o;))  = 


m(u;o  —  ol)2) 


y/m2(u;o  -  w2)2  +  (M2 

Tw 

x/m2(u;2-u;2)2  +  (ra;)2 


Tuj 


m(ui 0  —  u;2) 
one  finds 

m2(u;Q— a;2)2+(Ta;)2  «  4m2WQ((a;o-w)2+r2/4m2)  =  4m2u>o((l—  i4;/cvo)2+(T/2ma;o)2) 


If  o\  -&><§:  o\  one  finds 


If  we  define 


G{u>)  —  1/ -y/(l  —  lu  j  ujq)2  +  (T /2mu}Q)‘2 

we  see  that  if  we  take  the  square  of  it  we  get 
G(u>)  =  1/((1  —  w/w  o)2  +  (F/2mu;o)2)- 
So  the  width  of  the  lorentian  is 
_  _  f  _  _r  _ 

^ up  down 

m  Q 

and  we  get 


G{  uj)  —  1/ ^(1  —  +  (1/2Q)' 


tan(9?(o;))  — 


To; 

m(u o  —  u>2) 


1 

2Q(1  —  u/ujq) 


At  the  end  we  can  write 
z(t)  =  ^~G(u>)  cos(o )t  +  ip  —  y?(u/)) 


Since  we  have  excitation  from  the  source  and  from  the  drain  we  can  find  solution 
of  the  form  of  linear  superposition  of  the  two  independent  solutions,  meaning: 


f  f 

^(t)  —  G (u^)  cos(u;9t  T-  y>g  ^ ^ )  )  T-  ^ j—G(u>gd)  cos iyJsdf  sd 


The  expression  for  the  current  is: 

BC 

SI  =  7^r{SVg  COS(u Jgt  +  <Pg)~ 


dVf 

14,' 


z  ln(2z/r)  COS^9*  +  -  VffM)  +  cos(u;sdt  +  <A«f  -  y>«f(w))) 

414/  COs(uJ9dt  T  (psd)  ( '1 ) 

Let's  give  expression  for  fg  and  fS(j: 

0/7  ^  1/DC 

f  — _ lyDC§y  —  —C _ ^ _ £1/ 

•fa  V9  0V0  1  ' 0Vn 


fsd , 


dz 


,zln(2  z/r)  a 


BC  VDC 

fsd  =  ~^-VgDCSVsd  =  Cg-g-  sv* 
dz  y  z\n(2z/r) 

If  we  plug  these  expressions  and  use  the  following  relation,  cos(a)  cos(b)  = 

1/2  cos(a  +  b)  +  1/2  cos(a  -  b),  and  looks  at  the  low  frequencies  terms,  the  current 
looks  like 


1  OG 

4/  =  -  —SVg8Vsd{cos((u)g  -  usd)t  +  (pg-  ipsd)+ 

VDC  C 

(zln(2z/r)^  ^2fcG^  COS^  ~  “sd)t  +  *9  “  “  ¥>««) 

Cn 


2 •  G{uJsd')  COS^ipgd  ^sd{^)  V^sd)]} 


We  can  multiply  with  the  reference  signal  in  the  Lock-In  cos((cog  -coSd)t)  and 
perform  integral  over  few  periods.  It  is  important  to  understand  what  the  lockin  does: 
The  lockin  receives  the  reference  signal  and  we  always  have  the  freedom  to  say  that 
the  phase  is  zero.  Then  it  multiplies  the  input  signal  with  cos(cot)  and  sin(cot)  and 
perform  integral  over  few  periods  of  the  oscillations. 


If  we  assume  that  the  signal  is  A  cos(cot  +  qp),  one  finds: 


1_ 

T 

_1 

T 


fT  A 

/  A  cos(u )t  +  ip)  cos(u )t)  =  —  COs(<^) 
Jo  2 

fT  A 

/  A  cos(u>t  +  y?)  sin(u;t)  =  —  —  sin(y?) 
Jo  2 


Combine  these  results  together,  the  expression  for  the  current  will  be  given  by: 


1  r)C  VD( ”  C 

61  =  2dV95V9SVad{cOS{(p9  "  9sd)  +  (7H2^))22 k°M C°S(^  "  M 

If  Q  is  not  too  big  one  need  to  use  the  following  equations: 

GG(u>)  =  1/ \/(l  —  {uj/ujq)2)2  +  i^/^oQ)2 


tan((^(u;))  = 


Tuj 

m(ujQ  —  uj2) 


UJoUJ 

Q(uj2  -  uj2) 


1  dC  VDC  C 

61  =  2dVg6Vg6Vsd^C°S^9  ~  Vad>  +  ^z\n{2zlr)^~kGG^  COS^s  ~  ~ 

Having  these  formulas  one  can  fit  the  experimental  data  to  this  model.  The  result 
of  the  fitting  of  the  first  vibrational  mode  is  depicted  below: 


delta  phase  =  169.81 1 ,  Q  =  13.0157 


F[MHz] 


As  we  can  see,  the  agreement  between  theory  and  experiment  is  excellent.  The 
calculated  quality  factor  is  Q~13,  which  is  quite  low  for  this  setup.  Others  modes  have 


similar  Q,  up  to  30.  This  low  Q  is  a  result  of  the  suspension  process.  Currently  we 
developed  novel  approach  in  order  to  overcome  this  obstacle.  From  the  fitting  one  can 
extract  also  8z  and  the  applied  force.  It  turns  out  that  8z~10nm  and  8f~1  fN/VHz. 
These  values  are  not  the  minimum  detection  sensitivities  of  this  resonator  and 
measurement  of  these  values  are  still  in  progress. 


Atomistic  Simulations  of  Nano-Electro-Mechanical  Systems  Based  on  Carbon 
Nanotubes. 

The  primary  goal  of  the  proposed  research  is  to  take  advantage  of  recent 
technological  breakthroughs  in  the  fields  of  nanomechanics  and  nano-  and  molecular- 
electronics  in  order  to  study  suspended  nano-electro-mechanical  resonators  based  on 
Carbon  Nanotubes  (NTs)  [4],  We  plan  to  develop  a  theoretical  understanding  of 
dissipation  processes  and  nonlinear  dynamics  in  these  resonators  which  will  help  the 
experimentalists  utilize  them  as  ultra-sensitive  mass  and  force  sensors.  Specifically, 
we  intend  to  study  the  dependence  of  the  frequency  and  mode  shape  on  the  resonator 
quality  factor,  and  to  design  it  for  optimal  operation.  Currently,  the  conventional 
NEMS  have  the  advantage  of  having  higher  quality  factors,  higher  demonstrated 
resonance  frequencies  (  1GHz),  better  frequency  stability,  and  easier  on-chip 
integration,  while  NTs  possess  extremely-small  masses  and  therefore  hold  the 
potential  of  becoming  the  ultimate  mass  sensors  once  their  remaining  properties  have 
been  optimized. 

Results 

We  have  begun  atomistic  simulations  to  generate  data  for  the  SWCNT  vibrations 
in  order  to  find  and  tune  the  frequencies  to  obtain  improved  sensitivity  of  NT  devices. 
We  started  with  a  small  model  of  a  98.38  A  armchair  nanotube  with  a  diameter  of  9.4 
A  and  a  chiral  vector  of  (7,7)  clamped  at  both  ends  using  the  Brenner  interatomic 
potential  [5], 

Our  first  numerical  experiments  are  monitoring  thermal  vibrations  of  such  a  NT  at 
300K.  We  let  the  nanotube  vibrate  during  626  psec  with  a  time  step  of  0.5  fsec.  The 
center  of  mass  for  each  period  of  the  nanotube  (28  atoms  in  each  period),  see  Fig.  3, 
was  calculated  (red  dots/dark  gray  in  Figure  4),  we  call  this  line  of  points  the  Centered 
Axis  (CA).  We  analyzed  the  vibrations  of  the  CA  at  a  point  near  the  center  of  mass  of 
the  nanotube  with  clamped  ends  (3  periods  at  each  end  do  not  move).  The  tube  axis  is 
in  the  Y  direction,  X  and  Z  are  perpendicular  (see  Figure  4).  This  nano  tube  has  40 
periods  and  the  center  of  mass  (CM)  of  entire  tube  is  near  Y  =  20. 


one  period;  b) 


two  periods  of  a  (7,7)  Single  Walled  Carbon  Nanotube. 


Figure  4.  Single  walled 
Carbon  Nanotube 

(green/light  gray),  center  of  mass  (CM)  for  each  period  of  nanotube  (red/dark  gray), 
blue/dark  gray  atom  in  the  center  indicates  CM  of  the  nanotube. 

In  order  to  see  the  first  four  frequencies  (including  n  =  2  and  n  =  4)  we  monitored 
the  amplitude  near  the  CM  of  the  nanotube  vs.  time  of  the  vibration  (time  domain) 
and  from  the  Fourier  Transform  (FT)  analysis  the  amplitude  vs.  frequency  (frequency 
domain).  Longer  SWCNTs  were  tested  and  demonstrated  the  /  ~  L  2  dependency, 
where  /  is  the  frequency,  L  is  the  length  of  the  nanotube,  as  was  theoretically 
expected. 

First,  we  wanted  to  test  the  reliability  of  the  method  of  Fourier  Transform  for  the 
mode's  frequency  decomposition.  We  let  the  nanotube  vibrate  10  times  longer  than 
the  period  (T)  of  the  lowest  mode  (n  =  1).  The  vibration  of  the  CM  of  the  nanotube  in 
Z  direction  is  presented  in  Fig.  5  and  the  frequencies  in  Fig.  6,  this  figure 
demonstrates  that  the  frequencies  of  the  first  modes  can  be  distinguished  with  good 
resolution. 


Time  domain  of  Z  at  ¥  =20. 
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Figure.  5.  Time  domain  of  the  thermal  vibrations  of  the  CM  in  the  Z  direction  at 
300K. 

We  also  observe  in  Fig.6  that  the  amplitude  of  the  first  mode  (n  =  1)  has  a  high 
value,  however  the  amplitude  of  the  second  mode  ( n  =  2)  is  lower  than  that  of  the 


third  mode  (n= 3).  This  is  obvious  because,  since  the  center  of  mass  of  20th  period  of 
the  nanotube  is  near  the  CM  of  the  nanotube  (the  distance  being  about  1.42A),  which 
does  not  move  in  the  second  and  fourth  modes,  its  amplitude  is  low  (see  Fig.  7  for  the 


FFT  of  Z(t),  Y  =  20 


Frequency  (Hz) 

modes  of  vibrations). 

Figure  6.  Frequency  domain  of  the  thermal  vibrations  of  the  CM  in  Z  direction  at 
300K. 


n  =  2 


n  =  3 
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Figure  7. 

Vibration  and  PI 
standing 
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fundamental 

and  the  first  3  overtones  ( n  =  1,...,4). 


We  compared  our  results  to  frequencies  which  were  calculated  analytically  from  a 
continuum  model  [6]  as  presented  in  Table  1.  Since  we  treat  a  98.38A  nanotube, 


which  is  clamped  at  both  ends  (3  periods  at  each  end),  the  length  of  the  vibrating  part 
of  the  nanotube  is  98.38A  minus  length  of  the  6  periods,  which  do  not  move,  and 
equal  to  83. 62 A. 

Table  1.  Analytically  calculated  and  simulated  frequencies  in  THz.  t  is  the  wall 
thickness,  Eli  =  1.300  g/cm3  and  □fc  =9.517  g/cm3. 
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1 

2 

3 

4 

Simulated 

0.269 

0.659 

1.130 

1.660 

,  Pl 

0.279 

0.770 

1.510 

2.500 

£=0.34nm 

%error 

3.6 

14.4 

25.2 

33.6 

P2 

0.213 

0.586 

1.150 

1.900 

%error 

26.3 

12.5 

1.7 

12.6 

sidering  only  the  continuum  model  with  density  =  1.300  g/cm3,  which  is  a  more 
realistic  density,  the  first  mode  is  in  excellent  agreement  with  the  simulations.  The 
error  of  simulated  frequencies  vs.  analytically  calculated  ones  increases  with  n.  We 
still  need  to  investigate  this  effect,  but  think  that  the  reason  may  be  the  nanometric 
scale  effects  which  may  change  the  elastic  properties  of  the  SWCNT.  It  is  important 
to  note,  that  in  the  analytic  model  a  constant  value  of  the  Young's  modulus  of  1  TPa 
was  taken  for  the  calculations. 

To  be  more  specific  a  detailed  graphs  are  present  below. 

1.  We  have  calculated  the  vibrations  of  the  longer  (7,7)  nanotubes: 

R  =  4.72A;  nanotubes  of  83.62A,  125.43A,  167.24A  and  209. 05A,  40,  60, 


x10«  Dependance  of  Frequency  on  the  Length  of  the  SWCNT 


We  also  analyzed  the  error  (difference)  between  the  “old”  analytical  model 
and  my  results  and  plotted  a  graph  of  error  as  a  function  of  the  length  of  the 
nanotube  for  every  mode,  the  results  are  presented  in  the  figure  bellow: 


z(A)  Eror% 


Eror  as  a  function  of  the  Length  of  the  SWCNT 


The  vibration  of  Center  of  Mass  of  167.24A  nanotube  in  the  Z  direction: 

80p  microcan  255k:Time  domain  of  Z  at  Y=40 


2.  The  energy  changes  as: 
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3.  We  also  calculated  the  same  parameters  for  the  98.38A  nanotubes  with 
different  radii: 

(5,5)  R=  3.39A,  (7,7)  R  =  4.75A  ,  (9,9)  R  =  6.10A,  (1 1,1 1)  R  =  7.46A 
and  (14,14)  R  =  9.49A. 

The  results  are  presented  below: 


Future  Plans 

Our  next  step  in  the  future  will  be  a  test  of  different  boundary  conditions.  For 
example,  it  would  be  interesting  to  follow  the  frequency  changes  if  instead  of  the 
nanotube  being  clamped  at  both  ends  (3  periods  at  each  end)  only  the  atoms,  which 
are  supposed  to  be  connected  to  the  surface,  would  be  frozen.  In  later  stages,  in  order 
to  study  the  external  energy  loss  mechanisms,  we  will  connect  the  nanotube  to  the 


silicon  surface  at  the  edge  (see  Fig.  1  or  Fig.  8)  by  non-co valent  bonds  and  monitor 
the  changes  in  energy,  frequencies  and  the  quality  factor  of  the  system.  In  this  part  of 
the  research  electronic  effects  play  an  important  role  and  need  to  be  studied  using  ab 
initio  models.  Our  previous  experience  in  the  electronic  structure  calculations  could 
be  very  useful  for  this  purpose.  Some  additional  numerical  experiments  will  also  be 
carried  out.  To  study  the  influence  of  the  boundary  conditions  on  the  Q  we  will  follow 
the  changes  in  the  angle  between  the  nanotube  and  the  silicon  surface  as  is  shown  in 
Fig.  8. 


(a) 


Figure8.  Suspended  Single  Walled  Carbon  Nanotube,  showing  optional  angles  (a) 
-  (c)  between  the  nanotube  and  the  surface  (yellow  area). 

In  the  experimental  work  the  nanotube  was  suspended  between  two  electrodes,  the 
gate  voltage  was  applied  and  caused  an  electrical  force  downward  on  the  nanotube. 
Since  the  nanotube  has  an  extremely  high  Young's  modulus  ~  1  TPa  it  was  deformed 
as  a  result  of  the  applied  force  and  the  frequencies  of  the  vibrational  modes  were 
changed,  one  of  the  most  important  questions  is  how  the  frequencies  were  changed 
and  what  is  the  effect  of  such  a  change  on  the  quality  factor?  In  order  to  answer  these 
questions  we  will  study  the  tensioned  nanotube:  we  will  expose  the  nanotube  to  the 
external  forces  in  Z  and  Y  directions  as  shown  in  Fig.  9. 


(a)  (b) 


Figure  9.  Suspended  Single  Walled  Carbon  Nanotube:  (a)  external  forces  acts  in 
the  direction  opposite  to  the  Z  direction;  (b)  under  tension  in  Y  direction. 

Another  of  the  planned  actions  is  to  simulate  one  of  the  applications  of  the  Single 
Walled  Carbon  Nanotube  NEMS  by  adsorbing  a  small  molecule  onto  the  nanotube 
and  monitoring  the  frequency  shift  to  calculate  the  mass  of  the  attached  molecule. 


Suspended  Silicon  Nano  wires 

In  order  to  test  different  aspects  of  vibrational  modes  in  NTs  we  decided  to 
fabricate  also  silicon  nanowire  (SiNW)  resonators.  They  have  some  similarities  with 
NTs  and  can  be  tested  with  the  same  setups.  First  we  had  to  fabricate  such  field  effect 
transistor  devices.  We  study  how  to  make  optimal  metallic  contacts  to  these  wires, 
their  electrical  performance,  and  their  dependence  on  doping  parameters. 

We  have  studied  the  influence  of  the  RTA  procedure  on  our  pre-grown  SiNWs, 
finding  the  optimized  conditions  for  the  silicide  contacts  formation  in  a  systematic 
manner. 

Nickel  silicide  formation  is  one  of  the  most  effective  ways  to  control  and  improve 
the  contacts  in  silicon  NWs.  A  two  step  annealing  process  is  a  very  effective  way  to 
make  reliable  and  low  resistance. 

The  first  annealing  step  between  200-300  °C  forms  Ni2Si  with  low  diffusion  rate 
but  high  resistance,  while  the  second  step  between  400  and  450  °C  is  required  to  form 
low  resistance  NiSi  with  ~0.6  pm/sec  diffusion  rate  of  the  silicide  along  the  wire. 

To  explore  the  influence  of  the  RTA  temperature  on  the  transport  behavior,  the  two 
steps  RTA  approach  was  repeatedly  executed,  where  the  temperature  of  the  second 
RTA  step  was  raised  in  50  0  up  to  450  °C. 

To  estimate  the  impact  of  the  separated  RTA  steps,  Ron  (2  probes  channel 
resistance  taken  from  the  saturated  current  of  the  strongly  accumulated  NW)  was 
extracted. 

Typically,  for  a  constant  S/D  voltage  when  the  swing  voltage  of  the  gate  is 
considerably  increased,  the  current  start  to  saturate.  The  main  reason  for  current 
saturation  is  ascribed  to  the  serial  contacts  resistance  whereas  the  silicon  channel 
impedance  itself  is  relatively  negligible  for  such  a  high  gate  bias.  The  external 
contacts  resistance  can  be  estimated  by  plotting  K/s  I  Id  ~  Ron  ~  R  NW  +  Rext  as  a 
function  of  l/(Vg-VT).  Alternatively,  one  can  extracts  Ron  by  dividing  the  applied 
bias,  Vds,  with  the  saturation  current  which  is  acquired  at  the  highest  gate  voltage 
(figure  1,  red  circles). 

In  addition,  in  the  case  of  a  reliable  ohmic  contact,  one  can  expect  no  dependence 
between  Ron  and  Vds,  while  for  a  rectifying  contact  the  resistance  should  varied. 

Namely,  the  low  contact  resistance  should  maintain  for  a  large  range  of  Vds  biases. 

Under  these  criterions  more  than  35  devices  were  examined,  5  up-to  8  devices  for 
each  RTA  step.  Only  devices  with  clear  gate  dependence  that  can  be  firmly  opened 
and  closed  were  taken  into  the  statistics.  Figure  2  summarizes  the  results  of  the 
average  2-probes  resistance  versus  the  annealing  temperature  for  constant  bias  of 
Vds=  0.8  V. 
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Figure  1  :  Ids  Vs.  Vg  at  Vds=  -0.1  and  -0.5  V  for  different  RTA  procedures:  no  annealing,  200 
and  250  °C.  The  red  circles  mark  the  saturation  current  of  the  strongly  accumulated  NW,  for 
Ron  extraction. 


Figure  2  :  Ron  Versus  the  temperature  annealing  for  constant  bias  of  Vds=  0.8  V.  The  red 
bars  represent  the  standard  deviation  at  each  temperature. 


The  clearest  improvement  in  Ron  was  obtained  around  400-420  C°. 

Moreover,  the  results  exhibit  a  minor  variation  between  different  devices  with 
typical  resistance  around  250-350  KQ,  reaching  to  the  lowest  resistance  of  170  KQ. 


Further  increase  of  the  RTA  temperature  impaired  the  contact  performances  and 
deformed  the  NW  shape. 

Also,  we  discovered  that  for  a  short  period  of  RTA  (~10  to  15  sec)  at  420°C,  a 
reliable  NiSi  contact  was  formed  with  negligible  silicide  extensions  of  ~50nm,  which 
is  relatively  insignificant  in  compare  to  the  1  pm  NW's  channel  length. 

It  should  be  mentioned  that  the  interface  between  the  silicide  and  the  silicon  appear 
to  be  sharp  with  no  segmentation  or  nickel  contamination  in  the  SiNW  midsection. 

Accordingly,  15  sec  at  420°C  RTA  procedure  was  chosen  as  the  conventional 
annealing  process  for  manufactured  back-gate  devices. 

For  transport  measurements,  more  than  30  phosphorus-doped  (n-type)  SiNW 
devices  where  fabricated  and  their  electrical  characteristics  were  assessed.  Figure  3 
shows  typical  Ids  current  versus  drain-source  voltage  for  28  nm  diameter  and  1 .6  pm 
length  (in-between  the  S/D)  SiNW. 

Examination  of  individual  Ids  vs.  VdS  curves  shows  that  the  current  begins  to 
saturate  at  positive  but  not  negative  values  of  Yds  and  the  conductance  increases 
(decreases)  as  Vgs  becomes  more  positive  (negative). The  saturation  and  Vgs 
dependence  are  characteristics  of  accumulation-mode  n-channel  FET,  and  thus 
demonstrate  that  the  phosphorous  has  been  incorporated  as  active  donor-dopant  in  the 
SiNW. 

In  addition,  Ids  vs.  Vg  (figure  4)  curves  exhibit  high  on-currents  that  reach  to  more 
than  4  pA  and  relatively  high  transconductance  of  ca.  600  nS  (for  Vds  =  1  V). 
Furthermore,  the  SiNW  device  can  be  firmly  depleted  and  the  current  can  be  cutoff 
efficiently.  The  full  on/off  current  ratio  determined  from  logarithmic  plot  is  greater 
than  104  where  the  sub-threshold  slope  is  about  -1000  mV/decade. 


Figure  3  :  Typical  Ids  versus  Vds  curves  of  28  nm  diameter  and  1.6  pm  length  n-type  SiNWs, 
The  back-gate  voltage  for  each  curve  is  indicated. 


Figure  4  Typical  Ids  versus  Vg  curves  of  28  nm  diameter  and  1.6  jam  length  n-type  SiNWs, 
Vds  voltage  for  each  curve  is  indicated.  Insert  shows  the  logarithmic  scale  plot  of  Ids  versus  Vg 
measured  at  Vds=l  V,  presenting  on/off  ratio  >  104  and  sub-threshold  slope  of  -1000  mV/decade. 


Equivalently,  20  boron-doped  (p-type)  SiNW  devices  were  fabricated  and  their 
electrical  characteristics  were  assessed.  Typical  transfer  characteristics  are  presented 


in  figures  5-6. 
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Figure  5  :  Typical  Ids  versus  Vds  curves  of  36  nm  diameter  and  2.4  pm  length  p-type  SiNWs, 
The  back-gate  voltage  for  each  curve  is  indicated. 
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Figure  6  :  Typical  Ids  versus  Vg  curves  of  36  nm  diameter  and  2.4  pm  length  p-type  SiNWs, 
Vds  voltage  for  each  curve  is  indicated.  Insert  shows  the  logarithmic  scale  plot  of  Ids  versus  Vg 
measured  at  Vds=-1  V,  presenting  on/off  ratio  >  104  and  sub-threshold  slope  of  450  mV/decade. 


As  can  be  seen  in  figure  5,  the  current  begins  to  saturate  at  negative  but  not 
positive  values  of  Vds  and  the  conductance  increases  (decreases)  as  Vgs  becomes 
more  negative  (positive)  as  expected  from  the  p-channel  SiNW-  FET.  The  extracted  2 
probe  resistance  is  about  230KT2. 

The  Ids  vs.  Vg  curves  exhibit  high  on-currents  that  reach  to  ca.  5  pA  and  high 
trans-conductance  of  900nS  (for  -  IV  Vds  bias). 

The  subthreshold  slope  is  450  mV  per  decade,  which  is  a  great  improvement  in 
regard  to  the  slope  of  the  n-type  device.  As  before,  the  wire  can  be  firmly  depleted 
and  the  current  can  be  cutoff  with  no  observable  ambipolarity  behavior,  yielding  more 
than  104  on/off  current  ratio. 

In  general,  p-type  back-gate  devices  exhibited  better  consistency  in  their 
performance  and  provided  superior  results  than  the  n-type. 

To  realize  the  differences  between  the  two  NW-FET  types,  statistics  and 
histograms  of  the  NW's  key  parameters  were  taken,  where  green  and  red  colors 
indicate  the  n-type  and  the  p-type  frequent  results,  respectively.  In  addition,  two 
colored  arrows  in  each  figure  mark  the  average  empirical  result.  All  electrical 
characteristics  data  extracted  from  single  NWs  when  biased  at  |Vds|  =  1  V. 


Figure  7  SiNWs  diameters  Histogram,  extracted  by  AFM  measurements,  P-type  NWs  are 
mostly  in  between  30  to  40  nm  while  N-type  NWs  are  in  between  40  to  55  nm.  The  average 
thicknesses  are  34  and  48  nm  for  p-type  and  n-type,  respectively. 
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Figure  8  Histogram  of  Ni-contacted  SiNWs  resistance  determined  by  2-probe  measurements, 
the  average  resistances  are  430  KQ  and  640  KT2  for  p-type  and  n-type,  respectively. 


Figure  9  Histogram  of  SiNWs  mobility  valued  from  the  max  tarns-conductance  in  the  linear 
zone,  exhibiting  high  effective  mobility  with  average  of  91  and  65  cm2/V*sec  for  p-type  and 
n-type  NWs,  respectively. 


Figure  10  Histogram  of  the  transconductance  (gm)  taken  from  the  maximum  slope  of  the 
linear  zone  of  the  Id-Vg  curves.  The  average  results  are  820  and  475  nS  for  p-type  and  n-type 
NWs,  respectively. 


Figure  11  Histogram  of  the  max-current  density  (Js)  extracted  from  electrical  characteristics 
in  the  on-state  divided  by  the  NW  cross-section.  The  average  results  are  0.33  and  0.23  MA/cm2 
for  p-type  and  n-type  NWs,  respectively. 
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Figure  12  Histogram  of  the  sub-threshold  slope  for  thelOO  nm  oxide  back-gated  devices,  the 
average  slopes  are  545  and  1000  mV/decade  for  p-type  and  n-type  NWs,  respectively. 

The  superiority  of  p-type  devices  in  comparison  to  the  n-type  is  evident  in  all 
extracted  electrical  characteristics.  The  result  might  be  surprising  since  in  bulk  silicon 
the  electron  mobility  is  higher  than  hole.  But,  taking  into  account  the  NWs  diameters 
(where  p-type  found  to  be  ~15  nm  thinner  than  the  n-type  NWs  in  average)  and  the 
higher  built-in  Schottky  barrier  for  the  n-type  wires,  the  present  situation  can  be  more 
understandable.  What  is  more,  examination  of  the  growth  conditions  of  this  two  n/p 
NWs  batches  confirmed  that  the  doping  concentration  of  the  n-type  wires  are 


considerably  lower  than  the  p-type,  a  factor  that  can  cause  to  another  derogation  of 
the  electron  injection  ability  through  the  SB  contacts. 

To  summarize,  Table  1  presents  some  key  characteristics  of  our  devices  in 
compare  to  state-of-the-art,  back-gated  SiNWs  reported  in  recent  years. 


Table  1  :  Comparison  of  state-of-the-art,  back  gated  SiNW-FETs. 
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resistance 

The  electrical  characteristics  of  the  fabricated  back-gate  devices  are  quite 
resembled  to  the  previously  reported  achievements. 

The  new  SiNWs-FETs  exhibit  the  same  magnitude  of  high  on-current  and 
transconductance  which  are  essential  to  high-speed  electronics. 

In  addition,  transport  behavior  exhibits  high  free-charge  mobility  (~100  -  150 
cm2/V  sec)  which  is  comparable  to  the  technology  of  planer  SOI  (silicon  on  insulator) 
FETs. 

The  high  on-off  current  ratio  with  the  relatively  low  capacitance  of  the  NWs 

(~1  fF  per  1pm  SiNW)  suggested  that  these  NW-FETs  are  suitable  for  low-power 
applications. 

Equipped  with  SiNW  transistors  we  developed  novel  technique  to  fabricate 
suspended  devices  with  minimum  contamination  of  the  wires  themselves.  This 
approach  is  crucial  in  order  to  achieve  higher  quality  factors.  Below  are  pictures  of 
various  suspended  oscillators  that  we  have  fabricated.  One  can  identify  doubly 
clamped  wires,  bent  wires,  free  cantilevers,  etc. . . 


Tilted  SEM  images  of  different  crossbeams  and  cantilevers  made  of  SiNWs,  illustrating  the 
great  variety  of  structures  that  can  be  produced  by  the  new  fabrication  technique.  Different 
contacts  configuration  yields  diverse  built-in  strain  in  the  crossbeams. 

Having  these  resonators  in  hand  enable  us  to  study  influence  of  dissipation 
processes  that  are  mutual  to  NTs  and  SiNWs,  and  ones  that  are  of  different  origin.  We 
expect  that  boundary  conditions  and  anchoring  points  will  play  a  major  role  in  both 
resonators,  but  tension  and  phonons  distribution  will  have  completely  different 
influence. 

Static  Deformation: 

One  popular  method  to  characterize  Young  modulus  and  built  in  tension  is  to 
utilize  force  curve  measurements  with  Atomic  Force  Microscope  (AFM). 


Figure  1 :  Schematic  of  the  SiNW  beam,  clamped  by  an  embedding  metallic  contact,  subjected  to 
a  point  force  inflicted  by  an  AFM  tip.  A  The  beam  prior  to  bending;  B  The  beam  during  bending 
test.  The  total  descend  is  the  sum  of  the  NW  deflection  and  the  cantilever  deflection. 

Acquire  force  curve  measurements  yield  the  following  graph 
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from  which  one  can  extract  the  wire  Young  Modulus  and  built  in  tension.  For  the 
example  above,  the  residual  stress  is  lOOMPa  and  Y=170GPa,  in  good  agreement 
with  theoretical  values. 

Having  the  static  properties  of  these  wires,  one  can  estimate  their  dynamical  behavior, 
and  more  specifically  their  resonance  frequencies. 

Chemical  Functionalization  - 

An  important  ingredient  of  this  research  proposal  is  the  attachment  of  the  small 
mass  to  the  surface  of  NTs  or  SiNWs.  The  attachment  can  be  very  rigid  and  strong,  as 
covalent  bonding,  but  can  be  also  looser,  like  Van  der  Waals  interaction.  Covalent 
bonding  requires  intensive  research  even  before  the  final  device  is  available. 
Therefore  we  decided  to  use  amine  chemistry  in  order  to  attach  small  molecules  onto 
the  1-dimesional  wires.  NTs  contain  carboxylic  acid  (COOH)  groups  on  their 
surrounding.  Thus,  one  should  have  amine  group  connected  to  the  attached  molecule 
in  order  to  achieve  covalent  bond  between  it  to  the  NT.  For  that  purpose,  we  tested 
attachment  of  gold  colloid  to  the  surface  of  NT.  We  exploit  amine  chemistry  (utilized 
Dicyclohexylcarbodiimide  (DCC)  and  1-Hydroxybezotriazole  (HOBT))  and  attached 
gold  colloid  with  aminoethanthiol.  Figure  13  presents  the  chemical  flow  of  this 
reaction.  Since  the  gold  colloids  are  visible  with  SEM  and  AFM,  we  could  examine 
our  chemical  functionalization  with  great  details.  Figure  14  depicts  success  of  such 
chemical  attachment  over  the  exterior  surrounding  of  NTs.  Functionalization  of 
SiNWs  is  currently  under  process. 
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Figure  14:  AFM  images  of  gold  colloids 
(5nm  and  15nm)  covalently  attach  to  NTs. 
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